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By  J.  Pretsch 


With  the  aid  of  the  method  of  small  oscillations  which  was  used 
successfully  in  the  investigation  of  the  stability  of  laminar  velocity 
distributions  in  the  presence  of  two-dimensional  perturbations,  the 
excitation  of  the  unstable  perturbations  for  the  Hartree  velocity 
distributions  occurring  in  plane  boundary -layer  flow  for  decreasing  and 
increasing  pressure  is  calculated  as  a  supplement  to  a  former  report. 
The  results  of  this  investigation  are  to  make  a  contribution  toward 
calculation  of  the  transition  point  on  cylindrical  bodies. 
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I.  STATEMENT  OF  THE  PROBLEM 


If  one  wants  to  make  a  theoretical  calculation  of  the  profile  drag 
of  bodies  in  a  flow  for  a  certain  direction  of  air  flow,  one  must  know  - 
in  addition  to  the  pressure  distribution  -  the  position  of  the  transition 
zone  where  the  laminar  boundary  layer  becomes  turbulent.   The  separation 
point  of  the  laminar  layer  forms  a  rearward  limit  for  the  transition 
point  on  a  section  of  this  body  surrounded  by  the  flow.   It  lies  in  the 
region  of  the  pressure  increase  at  the  point  where  the  velocity  distri- 
bution in  the  boundary  layer  has  the  wall  shearing  stress  zero.   This 
separation  point  is  a  fixed  point  of  the  profile  in  the  flow,  the 
position  of  which  does  not  shift  due  to  a  variation  of  the  Reynolds 

number  Re  =  -^—     (JJm  =   velocity  of  air  flow,   t  =  chord  of  the  body)  . 

One  may  calculate  it  according  to  the  well-known  approximation  method 
of  Pohlhausen  which,  for  prescribed  pressure  distribution,  provides  for 
every  profile  point  a  velocity  distribution  of  the  boundary  layer  with 
a  certain  form  parameter  X      (X  =   -12  separation).  As  forward  limit  for 
the  transition  point,  one  may  take  the  stability  limit  of  the  laminar 
layer  with  respect  to  small  two-dimensional  perturbations  which  were 
calculated  in  references  2  and  3  according  to  the  method  developed  by 
W.  Tollmien  (ref.  1).   According  to  this  method,  there  exists  for  every 
form  of  velocity  distribution  in  the  boundary  layer  a  so-called  critical 

/u  S*\    / 
Reynolds  number  Re*  ...  =  -§ —     (Ua  =  local  potential  velocity, 

V  v  '^ 
\    •  cr 

5*  =  local  displacement  thickness)  below  which  all  perturbations  are 

damped;  its  value  increases  as  the  form  of  the  velocity  profile  becomes 

TT  FvJf- 

fuller  (with  increasing  X.)  .   Where  the  Reynolds  number  Re*  =  — - — 

formed  with  the  actual  potential  velocity  Ua  and  the  actual  displace- 
ment thickness  5*  exceeds  this  critical  Reynolds  number,  there  begins 
the  instability  of  the  boundary  layer.   In  contrast  to  the  laminar 
separation  point,  this  stability  limit  is  therefore  not  fixed  on  the 
profile  for  a  normal  pressure  distribution  but  travels  forward  toward 
the  stagnation  point  on  bodies  in  flows  with  increasing  Reynolds 
number  Re   (see  fig.  l). 

The  actual  transition  point  which  lies  between  these  two  limits  - 
stability  limit  and  separation  point  -  is  known  to  likewise  shift  forward 
with  increasing  Reynolds  number  Re.   That  it  lies  only  a  certain 
distance  behind  the  stability  limit  (as  shown  by  a  comparison  of  experi- 
mental transition  points  and  theoretical  stability  limits)  is  plausible, 
because  the  excitation  of  the  unstable  perturbations  starts  only  at  this 
limit  point  of  the  stability  and  must  obviously  have  attained  a  certain 
degree  before  the  instability  further  downstream  leads  to  the  breakdown 
of  the  laminar  flow  configuration. 
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For  this  reason  it  seemed  necessary  to  calculate,  or  at  least  to 
estimate,  for  the  velocity  distributions  in  the  laminar  boundary  layer 
in  the  entire  instability  range  of  the  latter,  the  excitation  of  the 
unstable  perturbations  as  well,  with  the  objective  of  making  an  improved 
calculation  of  the  transition  point  possible. 

For  the  velocity  distribution  on  the  flat  plate  in  longitudinal 
flow  (Blasius  profile),  H.  Schlichting  (ref.  k)    has  already  determined 
the  excitation  quantity  as  a  function  of  perturbation  frequency  and 
Reynolds  number  in  a  part  of  the  instability  range;  and  for  the  velocity 
distributions  in  the  region  of  pressure  increase,  W.  Tollmien  (ref.  5) 
has  explained  the  behavior  of  the  excitation  (in  first  asymptotic 
approximation  for  very  large  Reynolds  numbers)  in  a  very  general  manner, 
neglecting  the  effects  of  internal  friction. 


II.  THE  GENERAL  DIFFERENTIAL  EQUATION 
DESCRIBING  THE  PERTURBATION 


Since  the  bases  for  the  method  of  stability  investigation  have  been 
discussed  in  detail  in  an  earlier  report  (ref.  3)>  we  can  refer  to  the 
results  attained  there. 

Let  U(x,y)   and  V(x,y)  be  the  tangential  and  normal  components 
of  a  plane  steady  boundary-layer  flow;  let  x  denote  the  length  of  the 
arc  and  y  the  normal  to  the  profile  contour.   The  stream  function  of 
the  two-dimensional  perturbation  motion  which  we  superpose  on  this 
basic  flow  is  assumed  to  be 


■  /     .\     /    \  ia(x-ct)     /    n  Bit  i(ax-3r"t) 
i(x,y,t)  =  cp(x,y)e         =  rp(x,y)e  x  e 

c  =  cr  +  iCj_   Bj_  =  c^a   Br  =  Cra 


J 


(1) 


here  t  denotes  the  time,  a  the  spatial  circular  frequency  of  the 
perturbation,  the  real  part  cr  of  c  its  phase  velocity,  and  the 

imaginary  part  Cj_  a  measure  for  its  excitation  (cj_>  0)   or  damping 
(c-j_  <  0);  besides  Bi   is  the  logarithmic  increment  of  the  excitation 
of  the  perturbation  amplitude,  and  Br  the  circular  frequency  in  time 
of  the  perturbation. 
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If  we  substitute  the  motion  originating  by  superposition  of  the 
boundary -layer  flow  U,  V,  and  the  perturbation  motion 


u  =  **   =  q,'eia(x-ct) 


v  =  • 


Sy 


(2) 


5x 


iacp  +  &P\,ia(x-ct) 
dx; 


into  Navier -Stokes'  differential  equations,  we  obtain  -  as  was  proved 
in  reference  3  in  detail  -  the  differential  equation  describing  the 
perturbation  in  the  form 


(U  -  c)(V'  -  a2cp)  -  U"q>  =  — —  (q^V  -  2a2cp"  +  cAp) 
\         /  aRe* v  ' 

UaS* 


Re* 


J 


(3) 


Here  the  prime  ( ' )  denotes  differentiation  with  respect  to  the  wall 
distance  yj  the  velocities  are  referred  to  the  local  velocity  Ua  at 

the  boundary- layer  limit,  and  the  wall  distance  as  well  as  the  wave 


length  A  = 


2* 

a 


are  referred  to  the  local  displacement  thickness  5*, 


In  order  to  avoid  misunderstandings  it  should  again  be  emphasized 
that  -  in  spite  of  the  assumption  that  U,  V,  and  cp  be  functions  of 
the  arc  length  x  -  only  the  form  of  the  local  velocity  distribution  is 
decisive  for  the  stability  investigation  as  one  recognizes  from 
equation  (3).   The  immediate  effect  of  the  pressure  gradient,  however, 
is  negligibly  small  as  is  also  the  influence  of  the  x-dependence  of  the 
perturbation  amplitude  cp. 

The  boundary  conditions  of  the  differential  equation  (3)  result 
from  the  condition  that  the  perturbation  velocities  u,   v  vanish  at 
the  wall  and  that  the  friction  effect  at  the  outer  boundary -layer  limit 
(U' '  =  0)  has  disappeared.  With  their  aid,  the  calculation  of  the 
excitation  of  the  unstable  perturbations  may  be  reduced  to  an  eigenvalue 
problem,  the  solution  of  which  is  discussed  in  section  V. 

We  deal  first  with  establishing  the  particular  solutions  of  the 
differential  equation  (3),  limiting  ourselves  to  small  values  of  the 
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excitation  quantity  c-jj  thus,  the  general  solution  of  the  differential 
equation  describing  the  perturbation  (3)  may  be  represented  in  the  form 

v=l         v=l    V    U0    / 


There  cpv  denotes  the  particular  solutions  for  cj_  =  0  which  are 

obtained  in  the  calculation  of  the  limiting  curve  of  the  instability 
range  (neutral  curve)  in  the  a,   Re*-plane  (ref.  3)>  and  the  o>v 
signifies  additional  functions  for  c^  >  0. 

We  turn  first  to  calculating  the  integrals   cp-.*,   qv>*  and  the 
additional  functions  od-,  ,  CDp. 

III.  THE  SOLUTIONS  cpj*,  y  *      OF  THE  DIFFERENTIAL 

EQUATION  DESCRIBING  THE  FRICTIONLESS 
PERTURBATION  FOR  FINITE  EXCITATION 

If  aRe*   is  assumed  to  be  very  large,  the  differential  equation  (3) 
is  simplified  to  the  so-called  frictionless-perturbation  equation 

(U  -  c)(cp"  -  a2rp)  -  U"cp  =  0  (5) 


This  differential  equation  has  a  pole  of  the  first  order  at  the  location 
U  =  c  =  cr  +  ici  to  which  we  coordinate  the  point  yc*  of  the  complex 

y-plane.   In  the  neighborhood  of  this  singularity,  a  fundamental  system 
may  be  easily  indicated  by  series  development.   In  order  to  establish 
the  connection  with  the  case  of  the  purely  real  c  treated  before 
(see  ref.  3)>  we  first  give  the  relation  between  the  complex  yc*  and 
the  wall  distance  yc  of  the  "critical"  layer  U  =  cr. 

From 


u(yc)  =  cr  (6) 
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and 


u(yC*)  =  cr  +  ici 


there  follows 

u(yc*)  -  U(yc)  =  ici  =  (yc*  -  yc)u0'  +  .  .  .       (8) 
and,  with  limitation  to  the  terms  linear  in  c-^,  therefore 

y*  =  yc  +  ^L  (9) 

uo 

We  shall  now  indicate  the  construction  of  the  solutions  cp-i*,   rPo* 

of  equation  (5)  for  those  special  velocity  distributions  U(y)  by  which 
ve  shall  approximate  the  laminar  boundary  profiles  of  Hartree  (ref .  6) 
for  the  calculation  of  the  excitation  of  the  unstable  perturbations 
in  the  same  manner  as  we  did  before  in  reference  3  for  the  calcula- 
tion of  their  critical  Reynolds  numbers. 

(a)  Binomial  Velocity  Distribution 
(Pressure  Decrease) 

In  the  region  of  the  pressure  decrease,  we  used  the  approximation 
function 

U  =  1  -  (a  -  y)n     n  =  2,3A  .  .  .  (10) 

where  a  denotes  the  coordinate  of  the  point  of  junction  to  the  potential 
velocity. 

With  the  new  variable 

y  -  yc 


vl 


a  -  yc 


11 
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and  with 


cx-l  =  a(a  -  yc. 


(12) 


the  perturbation  equation  (5)  for  indifferent  perturbations  then  reads 


£ 


1  -  (1  -  yi)n 


d^P 


n-2r 


a^cp  +  n(n  -  l)(l  -  y^"^)  =  0 
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Using  the  relations  (9)  and  (11),  we  now  introduce  the  complex  variable 


yi*  =  : — —  =  yi  +  - — —  =  yi -* (ik) 


a  -  y. 


a  -  yc 


u0'(a  -  yc) 


With  the  abbreviation 


ci 


uo'(a  -  yc)     n(cr  -  i; 


(15: 


the  differential  equation  describing  the  frictionless  perturbation  for 
nondisappearing  excitation  then  is  transformed  into  the  form 


1  -  y,*)n  - 


nif^l  -  (1  -  y) 


n-1 


■ay 


2  *\ 


n(n  -  1, 


1  -  yi**n-2 


+  (n  -  2)  if -A  -  y{ 


*\n-3 


cp*  =  0     (16) 


in  which,  according  to  our  presupposition,  only  the  terms  linear  in 
f^  have  been  taken  into  consideration. 
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If  one  writes  equation  (l6)  after  multiplication  by 

n*2 


1.(1  -3r1*)I'  +  nlf1[l-(l-3r1*)1'-1] 


and  after  division  in  the   form 


d^tp 


y  *2  J2L  +  q>*  J_    PiV1   =  0  (IT) 


dy*2  !Zi 


the  first  solution  qp-j*  is  given  by 


yi*  52  ^i*"  =  — — K        L  <°l)  (l8) 


a   -  yc  V=0  a    -   yc\  uQ' 


with 


«>i  -  2_  §vyiv  (19) 

v=o 


The  series  coefficients  e  *  are  obtained  from  the  3V   according 
to  the  recursion  formula 


v=l 

v(v  +  l)ev*  +  ^  0   *e  *  =0    e0*  =  1  (20) 

u=o 


The  solution  cp2   for  all  y,*  is  found  to  be 


■Pi 


%*  -  V/^  *  - 1  CAW  *  -1*  r^ >»  '1*  <-' 
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In  this  equation,  the  logarithms  are,  however,  ambiguous;  one  must  cut 
the  complex  yy*-plane  in  such  a  manner  that 

3  *    JT 

-  i  «  <  arg  yi  <  _ 


If  we  put 


*v*yi*v  =  C  Wnv  +  ifi  YZ  *wiv 

V=2  V=2  V=l 


(22) 


cp2   may  he  represented  for  y-^  =   real  part  of  y-^*  >  0  at  first  in  the 
form 


^2*  =  1  +  XZ  Viv  +  ifi5Z  *wiv 


V=2 


V=l 


I1 "  "^  * lfi  £  wj 


In  y-L   + 


ifi 


lCj 


—    =  <P2 r  ^ 

7l)  u0 


(23) 


with 


av 


a  -  yc 


fcihyi   +i-^yl-lnyi)  + 


2ei in  yi  y~^y/ 

v=o 


(24) 
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In  the  transformation  of  equation  (21)  into  equation  (23)  the 
relation 


el*  =  el(X  -  ifl) 


(25) 


is  put  to  use;  it  is  easily  obtained  from  the  perturbation  equation  (5! 
if  U(y)   and  cp-^   for  c^  =  0,  and  for  c^  =  0  in  the  neighborhood 

of  the  singularity  U  =  c  are  developed  into  a  Taylor  series,  since 


2e?  =  (a  -  yc)V' 


V 


:26) 


and 


2en*  =  (a  -  yc). 


U0"  +  ici-° 


Un' 


U0'  +  iCl-2- 


V 


~  fa 


(a  -  yc) 


V 


Un' 


1  - 


icWUo" 


uo'\uo' 


=  (*  -  yc)-^—  (1  -  ifi] 

u0' 


(27: 


If  we  express  in  equation  (27)  the  exact  value 


V    U0'"s 


V   V 
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again  by  the  approximate  value  fj_,    this  step  is  motivated  by  the  fact 
that  due  to  the  small  value  of  c±,    one  deals  here  only  with  a  correction 


term  for  the  value  of 


Tl  '  ' 
^0 

U0' 


For  y,  <  0,   cpp*  is  represented  by  the  expression 


cp2*  =  1  +  £  b/y^  +  2e1*I*£-Ln 

V  =2  "  y°  \ 


yi 


+  1  ar 


g  yi 


cp2 


ICj 


f^ 


3        *  ^   « 

2  jr  <  arg  y*  <  - 


(28) 


with 


ov>  = 


a  -  yc 


2e1   In 


2e1cp]_/i 


2ijre, 


yi 


CPl 


V=0 


gvyl 


i*  -  yc  v=0 


Wlv)> 


(29) 


The  term  with  i  arg  y-,*   in  equation  (28)  was  obtained  by 

H.  Schlichting  (ref .  k)    and  W.  Tollmien  (ref .  5)  by  a  discussion  of  the 
general  perturbation  equation  (3)  in  the  neighborhood  of  the  singularity 
U  =  cr  +  icj_  in  a  similar  manner  as  the  "transition  substitution''  in 

the  critical  layer  for  purely  real  c  by  taking  the  friction  effect 
into  consideration. 


Since  we  are,  in  the  calculation  of  the  transition  substitution, 
concerned  with  a  representation  of  e-i   of  maximum  accuracy,  we  shall 

replace  e-i   by  the  expression  (26)  of  the  Taylor  development  of  the 

exact  Hartree  profile. 
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For  the  linear  and  the  parabolic  approximate  distribution 
[n  =  1,   n  =  2  in  equations  (10)  and  (l6)]  H.  Schlichting  (ref.  k) 
has  already  given  the  solutions  qv*,   qW*.   This  calculation  was  con- 
tinued for  n  =  3  and  n  =  h.      The  coefficients  3v*>  &v*  >     W* > 
s,      hy   are  compiled  numerically  in  table  1.   For  the  convergence  of 
the  power  series,  the  reflections  made  in  reference  3  are  valid. 

(b)  Sinusoidal  Velocity  Distribution 
(Pressure  Increase) 

In  the  region  of  pressure  increase,  the  Hartree  velocity  distri- 
butions in  the  boundary  layer  (see  ref.  3)  were  approximated  by 
formulation  in  terms  of  a  sine  formula  introduced  by  W .  Tollmien 
(ref.  5). 


U  =US  +  (1  -  Us)sin[?-=-||]  (30) 


where  s  denotes  the  wall  distance  of  the  inflection  point. 
With  the  new  real  variable 


IL^z  (31) 


and  with 


y2   2  a-  s 


CC2  =  -(a  -  s)a  (32] 


the  perturbation  equation  (5)  for  the  indifferent  perturbations  then 
read 


sin  y2  -  y2s  +  sin  y2s 


d 


2, 


2., 


—  -  o^rpj    +  sin  |y2  -  y2s)9  =  0  (33) 
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where  y2   was  put  equal  to 


?2i 


n   s  -  yc 


2  a  -  s 


(3*0 


We  now  introduce  the  complex  variable 


*2 


*  _  «  y  -  yc 


2   a-s 


y2   +  o 


*  yc  -  yc 


2   a-s 


=  y2 


ic- 


Un'  2(a  -  s) 


!35) 


With  the  abbreviation 


U0'  2(a  -  s)      (1  -  Us)cos  y2£ 


(36) 


the  frictionless  differential  equation,  describing  the  perturbation 
for  non-vanishing  excitation,  then  is  transformed  into 


sin 


v*yp*2 


-  y2s 

a22cp' 


s  ,   -   cos  y2s 


d29* 


+   sin  y2g   -   if^cos  fy2      -  y2 

(y2*  -  y2s)-if2  cos(y2*  -  y2s) 


sin 


rp*   =  0 


(37) 


If   (37)    is  written  like  equation   (l6)    in  the   form 


y2*2     *^L  +  y*f:   PiW1    =   0 


dy2 


*2 


i  =  l 


(38) 


the  solutions  rp-j.  j   <P2*  are  determined  by  analogy  with  the  equations 
(18),  (23),  and  (28).   One  has 


cpi   =  y2  z_  ev  y2 

2(a   -    s)  v=0 


IC-i 


CPl 


2(a   -    s)lri        U0 


7^1 


(39) 


11+ 
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with 


CO 


i =  2ZZ  * 


v=0 


VJ2 


(to) 


For     y?   >  0,    the   solution     cp2      is   obtai 


ned 


q>2*   =   1   +  )        bv*y2*v   +  2e^ 
v=2 


ircp^ 


2(a   -   s) 


In 


y27 


v  <  arg  y  *   <  _ 


(hi) 


If  one  puts 


*ir_*V     =    \  K     „    V 

v=2 


bv^2      +   if2  Z_     ^2 

v=l 


(to) 


cpo*      for     y2   >  0     may  be  represented   in  the   form 


00  00  . 

cp2*    =   1    +   Yl  V2V    +   if2  XI  V2V    +   ^lf1   + 
v=2  v=l  V 


if2  \ 


sin  y2s   cos  y2s 


jrcp1 


2(a   -   s) 


+    i*2   2Z  Svy2V)(ln  y2    + 


V=0 


if2' 


lci 
Po «>p 

V 


(to) 


with 


C0o    =    1 

2        2(a   -   s) 


i    h^     +      2fleiCPl      /JL   +  ln   ^ 


V=l 


2(a   -   s)  \y2        sin  y2s   cos  y2< 


2ei  ln  y2  H  svy2' 


(to) 


v=0 
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In  transforming  (kl)    into  (U3),  one  makes  use  of  the  relation 


e-L*  =  eifl  + 


Ifr 


sin  y2s  cos  y2s 


(V>) 


which  follows  from  the  two  equations 


.  a  -  b.Uq1 


JT     U, 


0 


(1*6) 


and 


'a  -  s)  V  +  ici 


ei' 


V" 

U0' 


UQ-  +  iCi 


uo" 
V 


(a  -  s)  U0" 
n  Uo' 


t  t        TT    T  I  ! 


1  - 


iCi/U0"   u0 


Uo'\u0'   u0' 


(a  -  s)  U0"  /  if 


tt     Ur 


1  + 


sin  y2s  cos  y2 


(hi) 


Cfu*   for  y2  <  0,  is  represented  by  the  expression 


cp2*  =  1  +  Z_  bv*y2 
v=2 


*„  *v     * 


(a  -  s) 


9  *  In 


y2   +  i  arg  yg 


92 


ICh 


Ur 


OV 


I  n  <  arg  y2*  <  2  \(kQ] 
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with 


>     av 


2  a  -  s 


$Ti    ^  a  -  sly2   sin  y2s  cos  y2s  I 


2el  ln|y2  Z_  6vy2 
"  v=0 


Svy2 


v=0 


+  2irreiJ 


Jtrp1 


2(a  -  s)  sin  y2£ 


cos 


*2t 


(*9) 


For  the  calculation  of  the  transition  substitution,  we  shall 
replace  the  term  ej_,  as  in  section  III  (a),  by  the  accurate  value 

(equation  (k6))    of  the  Taylor  development  of  the  exact  Hartree  profile. 

In  a  comparison  of  the  relations  (2U)  and  (hh)    or  (27)  and  (h-'j)    or 
(29)  and  (U9),  it  is  striking  that,  in  the  expressions  for  e^*  and  o^ 

for  the  sinusoidal  boundary -layer  profiles  in  case  of  pressure  increase, 
the  product  sin  y2s  cos  y2   appears  in  the  denominator  of  several 

terms.   The  sign  of  this  product  is  negative  when  the  inflection  point 
of  the  velocity  distribution  lies  more  closely  to  the  wall  than  the 
critical  layer  (V2s  <  0);  the  critical  layer  thus  lies  in  the  part  of 

the  velocity  profile  showing  concave  curvature.   The  sign  of  the  expres- 
sions divided  by  sin  yQ  cos  yQ   then  is  the  same  as  the  sign  of  the 

corresponding  terms  in  the  solutions  for  the  velocity  distributions  in 
case  of  pressure  decrease  which,  as  is  known,  have  concave  curvature  at 
every  wall  distance.   If,  in  contrast,  the  critical  layer  is  located 
between  point  of  inflection  and  wall  (V2s  >  0 j,  then,  in  the  part  of  the 

velocity  profile  having  convex  curvature  the  signs  discussed  are  reversed. 
In  that  case,  when  the  critical  layer  shifts  to  the  point  of  inflection 
itself  /Vps  =  0\ ,  the  behavior  of  e^*     and  cop   is  regular  since  then 


ic- 


ifr 


2  1 


Us) 


(50) 
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and 


03 


2(yC  =  s)   2(a  -  s) 


v=l 


hvy2A 


2(a  -  s 


i— -  In  y2 


(* 


>  0 


(51) 


u>, 


2(yC  =  s)  "  2(a  -  s) 


JS. 


V=l 


V 
V^2   + 


"Pi 


2(a  -  s) 


In 


y2 


.  ? 


2  a  -  s 


1 


y2  <  o^  (52) 


The  series  coefficients  Bv*,   ev*,   bv* ,   gv,   hv  for  the 

sinusoidal  basic  velocity  are  given  numerically  in  table  2.   For  the 
convergence  of  the  series  developments,  the  explanations  of  the  earlier 
report  (ref.  3)  are  again  valid. 


IV.  THE  FRICTION  SOLUTION  cp^*  FOR 
FINITE  EXCITATION 


Besides  playing  a  role  in  the  critical  layer  U  =  c,  the  friction 
of  the  fluid  is  of  importance  also  in  the  neighborhood  of  the  wall  where 
it  occasions  two  more  solutions  cp^*,   cph*  of  the  general  differential 

equation  describing  the  perturbation  (3)- 
If  we  introduce  the  variable 


y  -  yc 


(53) 


with 


e*  = 


[U0'  +  ic- 


1  u0- 


•1/3 


(5M 
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we  obtain  from  equation  (3)  in  the  limiting  process  e  — »0  for  large 
Reynolds  numbers  the  differential  equation 


iaS> 


dT) 


dT] 


(55: 


The  solutions  of  this  differential  equation  are  just  as  independent  of 
the  form  of  the  velocity  distribution  U(y)   as  was  the  case  for  the 
solutions  cpo,   cpii  for  the  excitation  zero  (ref.  1) ;   they  read 

corresponding  to  reference  2: 


V 


*3,V 


dV 


^  * 


1/2 
T]    H 


(D,(2) 


1/3 


3V">  J 


\3/2 


dTi 


(56) 


where  H   '     signifies  the  Eankel  function  of  the  first  and  second 
kind,  respectively. 

Since  the  Hankel  function  of  the  first  kind  increases  for  large 
wall  distance  beyond  all  limits,  it  cannot  be  contained  in  the  general 
integral  (k)    so  that  we  there  may  put  C^  =  0. 


V.  STATEMENT  AND  SOLUTION  OF 
THE  EIGENVALUE  PROBLEM 


After  having  found  the  particular  solutions  of  the  differential 
equation  describing  the  perturbation  (3)>  we  now  state  the  eigenvalue 
problem,  (which  results  from  the  boundary  conditions  of  the  perturbation 
equation)  for  investigation  of  the  excitation  of  the  unstable  pertur- 
bations. 

At  the  wall,  the  tangential  as  well  as  the  normal  component  of  the 
perturbation  velocity  disappear;  thus  one  has 


Cl*^!**  +  c2*cP2*w  +  c3*<P3*w  =  ° 


(57) 


clW  +  c2*rP2V  +  c3*?3*v'  =  0 


(58) 
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At  the  point  of  juction  y  =  a  to  the  region  of  constant  velocity, 

U"   is  zero  and  therewith  cp*"  -  a2<p*  =  0,  that  is  cp*  =  e"ay,  so  that 
the  third  boundary  condition  reads 


cl**l*a  +  c2**2*a  =  0 


(59) 


with 


V  a   rv   a 


*»'  +  °Va'   V  =  1>   2> 


(60) 


A  term  *o  *  does  not  appear  because  the  particular  friction  solution 
cpo*  has  already  been  damped  at  the  point  of  junction. 

In  order  that  the  three  homogeneous  equations  (57)?  (58),  and  (59) 
may  have  a  solution  different  from  zero,  the  determinant  must  disappear 


¥w     ^2  w     <P3  w 
<PlV    <P&V    <P3*v' 


*l*a     *2*a       ° 


=  0 


(6i: 


The  solution  of  this  determinant  yields  the  equation 


^V    <feV*l*a  "  *lV*2*£ 


(62) 


or  with  Schlichting' s  abbreviations 


1    93V 

e  no  ^3  w 


n  * 
^0 


=  F^q*)  =  F* 


(63) 
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<P2*v*l*a  "  <PlV2ja 


E*(a,cr,ci)  (6k) 


Vc*   <P*v'®l*a  "  ?l*w  V 
finally 

F*  =  E*(a,cr,ci)  (65) 

Therein 

no    =  -  —  (66) 

The  complex  equation  (65)  is  equivalent  to  two  real  equations  in  which 
the  parameters  a,   Re*,   cr,   and  c±     are  contained.   If  one  limits 

oneself  to  the  case  of  indifferent  perturbations  (cj_  =  0),  one  obtains 
from  them  after  elimination  of  cr  a  relation  between  the  Reynolds 

number  Re   and  the  wave  length  —  of  the  perturbation.   That  is  the 

equation  of  the  neutral  curve  or  the  curve  of  constant  zero  excitation 
by  which  the  stable  and  unstable  perturbation  states  are  separated  and 
the  lowest  Reynolds  number  of  which  is  the  so-called  critical  Reynolds 
number  Re*cr;  when  this  critical  Reynolds  number  is  exceeded,  the 
excitation  of  the  unstable  perturbations  begins. 

These  neutral  curves  were  determined  in  a  previous  report  (ref .  3) 
for  the  velocity  profiles  for  decreasing  and  for  increasing  pressure 
calculated  by  Hartree. 

In  the  present  investigation,  we  shall  assume  the  excitation 
quantity  c^  to  be  different  from  zero  and  investigate  the  curves  of 
constant  excitation  enclosed  by  this  neutral  curve. 

For  this  purpose,  we  shall  solve  the  complex  equation  (65)  in  a 
somewhat  simpler  manner  than  Schlichting  (ref.  k) ,  by  determining  not 

B  C-!  3  C,' 

the  differential  quotients and  —  at  the  location  of  the 

d  a       d  Re* 
neutral  curve,  but  the  curves  of  constant  excitation  cj  >  0  directly 

by  a  method  similar  to  that  applied  for  the  neutral  curve  c^  =  0; 

however,  without  making  use  of  this  curve  itself.  We  consider  first 
the  left  side  of  equation  (65). 
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Since  we  intend  to  limit  ourselves  to  small  values  of  excitation, 
we  develop  F*   in  the  form 


'(v)  =  fM  +  (^  - 


^0 


dF 


VV/V^o 


(67) 


where     t\~     as  before   in  reference   1   is   defined  by  the  relation 

,\  +1/3 


tio   =   -  yc(oRe*U0'^ 


(68) 


According  to  equations  (9),  (5*0,  (66),  and  (68)  one  now  has 

1/3 
10* 


c\       u0'yC; 


o^^Yi  +  ffiSalT 


-  ycfaRe*U0lN)l/3 


^0 


V2 


ici  /-,  L  yc  u0' 


1    +   i-  1    + 


U0'yc\         3    U0' 


1+^L/1  +  ^uo"' 


u0'yc 


3    uo' 


Thus  equation  (67)  becomes 


^^~*M^o(-r^ 


(69) 


=  G(Tl0'cr'ci) 


(70) 
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One  sees  therefore  that  the  form  of  velocity  distribution  which  does 
not  enter  into  the  exact  solution  F(t}0*  j  according  to  equations  (63) 

and  (56)  does  appear  explicitly  in  the  development  of  F^tjq*  \  with 

respect  to  c^. 

The  differential  quotient 


dF 


dF 


dFi 
+  1- 


&T\0  d.T\Q  &T\Q 


was  determined  by  graphic  differentiation  of  the  function  f/t]0).   The 

numerical  values  of  its  real  and  imaginary  part  are  given  numerically 
in  table  3«   Since  real  and  imaginary  parts  nowhere  disappear  simultane- 
ously, the  function  Ff^o*)   will  be  free  from  singularities  in  region 

around  the  function  F(t)q)  and  the  development  (67)  will  be  thus  per- 
missible. 

We  now  consider  the  right  side  of  equation  (65)  which  is  defined 
by  the  relation  (6k). 

Remembering  the  splitting  of  the  frictionless  solutions 
<Pv*(v  =  1,2)   in  the  form  of  equation  (k) ,   we  first  put 


with 


furthermore 


*  icj_ 

V   a   =  *va   "  ~7~T  ^Va 
u0 


^va  =  %a      +  c 


Va 


ICi 


(TD 


Vv  w  =  ^vw  -  — r  «\<w 


KV    W 


=  9- 


vw 


IC-i 


Ur 


OJ 


VW 
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For  the  binomial  velocity  distribution  U  =  1  -  (a  -  y)n  the 

values  of  the  additional  functions  co...   co,,   and  their  derivatives 

va        vw 

cov      ,     aov'        are  given  by  the  equations 
a  w 


co 


CO 

v=0 


CO, 


'2a 


a  -  yr 


v=i  a  -  yc 


CD 


la 


£ 


a  -  yc  v=i 


Vg, 


CO, 


'2a 


vhvi  + 


2ei 


2  9l£ 


(a   -  yc)2  v=1  (a   -  yc)2\       a   -  yc 


+  cpia     +  col£ 


COn 


iw  =  2^  Sv^ 


v=0 


CO, 


2w 


a  -  yc 


2e      In 


V=l 


lw 


V 


v    L  Sei^iw/  1 


lw 


a  -  yc  Uiw 


-  ln|yiw|^iw\  + 


CO 


lw 


+   ifl  ^0_/_Plw_ 

u0'  \a  -  yc 


-    CO 


lw 


>(72) 


CO 


lw 


a   -  y 


00 


c  v=l 


CO, 


2w 


00 

y-K  v=I 


2ei 


(a   -  yc)2   v 

,/   1_ 

lw 


(a  -  yc)2 


<Plw 


'iviy 


In 


yiw 


CD 


lw 


+  _Ji  +  in 

yiw 


yiw 


a  -  yc\     yiw" 

colw'(a  -  yc; 


yiw/ 


I  I   /  „       ' 


i)rHQ_/_flw_ 

U0'  la   -  yc 


co 


lw 


2k 
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For  the  sinusoidal  velocity  distribution 

U  =  Us  +  (1  -  Us)sin(y  '  s  *' 

la  -  s  1 

in  contrast,  the  corresponding  formulas  read 


*la  =  *—  ^2 
V=o 


u2a 


2(a   -    s) 


H   V2. 


v  +   "ei^Pla  /  1 


v=l 


In  y2s 


2el   ln  y2amla 


ula 


2(a   -    s)    v=l 

2  DO 


vgvy2£ 


v-l 


v-l  + 


if^ei 


=  MlTT^  fe  vhvy2a-     -  ^^r^ 


"qpia 


2(a-s)l      y2a2 


y2a  sin  y2s  cos  y2s 


sin  yo=   cos 


y2E 


y2s   cos  y2. 


">la  2(a   -   s)        -, 

—  „  +  ln  ^a^la 

y2a  * 


£ 


ulw  =  Z_  S-/2w 
V=0 


2w        2(a   -   s) 

2ei  ln  |y2w|a>i 


t "-"  *  ^fe  ♦ = 


In     y2w 


U0' 


y2s  cos  y2s 
[<plw 


Uri'\2(a   -   s)sin  y2s   cos  y2E 


^lw 


ulw 


u2w 


2(a   "   s'    v=0 

,2 


X!  vgvy2w 


v-l 


Ma   -   s)2^TX 


Hvhvy^-1.-^! 


fl2e 


2(a  -   s)< 


"rPlw 


2(a   -   s)  \       y2w: 


y2u   sin  y2g   cos  y2s 
^lw        2(a   -   s) 


+  »,'(+-+      ln  M \  + 


y2w 


In  y2wUlw' 


y2w       sin  y2s   cos  y2g/ 
._  U0"/  ■   "9lw 


f  (73) 


U0  '  l2(a   -   s)sin  y2s   cos  y2E 


">lw 
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If  the   critical  point   shifts   to  the   inflection  point  of  the  velocity 
distribution      (yc   =   s,   y2S   =  CM,    the  equations   for     u>>a>     ^2&   >     °^2\i> 

^2w'       (equation   (7*0)    are   simplified  to 


CD, 


2a 


2(a  -   s) 


ncpl£ 


CD, 


ff2 


2a 


k(a  - 


YL  hvy2aV    +  , 

V=l  2(a   -    s. 


'VoaV_1   + 


in  y2a 


jrcp1 


a 


V=l 


8U    •M^  +  ^     ln  ^ 


U3, 


'  2(a   -    s) 
J2. 


CD, 


2v 


Ma   -   s)2 

i7r3?lw' 
Ma   -   s)2 


5  W  *  sit^t  ta  ^ 


V2-   Vh  y      V-X   +  ^Iw 

fci    V2w  2(«  "  3)y2w 


1  tt3 


tPlw 


>(lh) 


Ma   -   s)2 

+  <Plv'  ln  |y2w| 


We   now  split     cp2w,      cp2w'  ,     °>>w>     ^w'      into  real  and  imaginary 
parts   according  to  the   formulas 


cp2w   =  A1    +   iBL  032w   =  M1    +    iNx 


P2w     =  A2   +   iB2  co2w     =  M2    +   iN2 


(75) 


then  one  obtains 


^  =  Ai  +  U7Ni  +  iBi-U7Mi, 


> 


qfeV    =  A2   +  ^  N2   +   i   B2    "  ^r  M2 


J 


(761 
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Then  we  put  the  expression  (6k)    into  the  form 


E*(a,cr,ci)  =  E(a,cr) 


(zl  +  iz2]j 


1  +  c^z-^  +  iz 


(77! 


where     E,    as  previously  defined  in  reference   1,    is   defined  by 


E   =    - 


1       ^w^la   "   cPlw*2a 
yc  <P2v'*la   "   ^lw'^c 


(78) 


First,  limiting  oneself  to. the  terms  linear  in  Cj_  one  has 


<P2Vl*a  ~  Via   +  ^(Nl°la  +  Via) 


Via   "  ^(Vla  +  Ml*la 


cpl*w*2*a  ~  <plw*2a   "  fj-r(m1^2a   +   cPi^a 


^vVa  -  Via  +  ^Vla   +  ^la)  + 


B2*la   "  -^r(A2%a   +  Via) 


Vw'Va  ~  rPlw'$2a   "  ^(^lw'*2a   +  ^lw'^2 


>        (79) 
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Hence,  there  results  after  a  short  intermediate  calculation 


E*   ~  E 


1    - 


ic< 


u0' 


u0  yc 


(Ni$ia   +  Bjflia)    +    i(cuiw<&2a   +  rPlw^2a   -  Ajflia    ~  Mi*ia) 


Ai*la    "   cPiw$2a   +   i*l9lt 


(Ng1>ia    +  B2nla)    +   iftlw'^a    +   ?lw'"2a    -  A2^la    -  M2^l£ 


A2$la   "  (Plw'*2a    +   iB2$la 


(80)' 


We  put  for  abbreviation 


Nl°la  +  Bl%a  =  ml 


a)lw*2a  +  ^lw^a  -  Alfila  -  Mi*ia  =  ni 


N2*la   +   B2%a   =  m2 


03 


lv'*2a   +  "Plv'^a   "  A2%a   -  M2*ia   =   n2 


Via   "  fPlw°2a  =  Kl 


A~<J> 


2°la   "   <Plw  °2a   =  *2 


"1 


> 


(81) 
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Then  there  follows  from  equation  (80) 

m-L    +   in-[_ 


E*~    E 


=  E 


1    - 


VYc        Uo' 


0 


m2    +   m2 


L 


Kl    +   iBl°la       K2   +   iB2*la 


1C 


^   +    ^r 


Uq  yc     u0 


mlKl   +   nlBl*la        n^Kg    +   n2B2$]_, 


KX2   +  B!2$la2  k22    +  B223>la2 


Id 


UO 


n]_K]_   -  m^B^ia        n2K2   -  mjjB^ia 


Ki2  +  Bi2$ia2    K22  + 


B22$ia£ 


(82) 


If  we  finally  introduce   the   designations 


Ur 


ml^l   +   nlBl*la       m2K2   +   n2B2(I)]_a 


Kl2    +   B^ct^2 


K22    +  B22d,la2 


^-  +   ^ 


u0  yc 


u, 


nlKl    -  mlBl*la        "2*2    ~  ^2°  la 


Kl2   +  B^Oia2  K22    +  Bp^la2 


(83) 


the  representation  (77)  is  attained. 

By  means  of  the  equations  (70),  (11),    and  (83)  we  thus  have  divided 
the  two  functions  F*  and  E*  (a,cr,ci)  appearing  in  equation  (65)  into 
real  and  imaginary  parts  and  can  now  graphically  solve  this  equation. 
For  this  purpose,  we  plot  for  constant  ci  for  several  values  yc  the 
imaginary  part  of  the  function  F* 


l(F*)    =1   F^o) 


1)0 


u 


1  +  ^VW 


0  yc 


3        U0'/dT)0 


(8U) 


against   its  real  part 


b(f*)  =b[f(^0)]   -no 


+  yc  U0      Wi 

uo'yd        3    Uo'idno 


(85) 
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and  plot  in  this  polar  diagram  to  the  same  values  c^  and  yc,  the 
imaginary  part  of  the  function  E*(co,cr,ci) 

ife*)  =  I(E)(1  +  ciZl)  +  ciz2R(E)  (86) 

against  its  real  part 

r(e*)  =  R(E)(1  +  ciZl)  -  ciz2l(E)  (87) 

The  points  of  intersection  of  the  functions  F*  and  E*  connected  one 
to  the  other  by  the  same  pairs  of  values  cr,   Ci  yield  first  the  values 
t]q  pertaining  to  the  corresponding  a-values,  and  then  with  (68)  the 
curve  a(Re*)   of  the  constant  excitation  c^.   This  plotting  and  cal- 
culation then  are  repeated  for  other  values   c^ . 

VI.  RESULTS  OF  THE  CALCULATION 


As  immediate  result  of  the  rather  extensive  calculations,  the  polar 
diagrams  with  the  curves  E*(a,ci,cr)   and  F*  =  F(t|q*)   for  several 

Hartree  velocity  distributions  U(y)   in  the  laminar  boundary  layer  are 
represented  in  figures  2  to  21.   These  boundary -layer  profiles  belong 
to  the  special  class  of  pressure  distributions 

Ua  =  const.  xm  (88) 

and  are  rigorous  solutions  of  Prandtl's  boundary-layer  equations  which 
were  obtained  with  the  aid  of  a  Bush  apparatus  (ref.  6).   We  use, 
according  to  Hartree,  for  characterization  of  these  boundary-layer  pro- 
files the  parameter 

2m  (89) 


1  +  m 


The  profiles  with  m  >  0,   3  >  0  occur  in  case  of  decreasing  pressure; 
the  profiles  with  m  <  0,   3  <  0  in  case  of  increasing  pressure.   The 
profile  3=1  is  the  exact  Hiemenz  profile  at  the  stagnation  point  of 
a  cylindrical  body,  the  profile  3=0  is  the  exact  Blasius  profile  at 
the  flat  plate  in  longitudinal  flow,  and  the  profile   3  =  -0.198   is  a 
separation  profile  (wall  shearing  stress  zero). 
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is  the  dimensionless  wall  distance  used  by  Hartree  and  the  coordinate 
a   of  the  point  of  junction  to  the  potential  velocity  Ua   is  connected 
with  the  displacement  thickness  by  the  relation 


a_   =  £aV|2  -  3 
5       k6 


(92) 


the  quantity  kg*(3)  may  be  found  plotted  graphically  in  figure  2  of 
reference  3- 

The  wall  distance   £c   of  the  critical  layer  indicated  in  figures  2 
to  21  is  obtained  if  in  equation  (91)   y   is  replaced  by  yc . 

The  following  basic  remarks  should  be  made  concerning  the  graphic 
solution  of  the  eigenvalue  problem  in  figures  2  to  21: 

Since   c-j   has  been  presupposed  so  small  that  only  the  variations 

of  the  particular  solutions  cpv*  linear  in  c-?   need  to  be  taken  into 

consideration  compared  to  the  solutions  cpv  for  c^  =0,  the  curves  E* 

and  F   appear  for  equal  c^-interval  as  "equidistant"  curve  families 

in  the  sense  of  equations  (77)  and  (70).  Actually,  this  "equidistance" 
will  be  lost  in  case  of  higher  values  of  c^;  however,  the  calculation 

expenditure  would  increase  intolerably  even  if  only,  for  instance,  the 
terms  quadratic  in  c±      were  to  be  taken  into  consideration.   In  this 
sense,  the  curves 

Pi 8*   c±ab* 

~vT  =  ua 

against  a&   represented  in  figures  22  to  26  (from  now  on  we  use  dimen- 
sional quantities)  which  result  from  the  evaluation  of  the  polar  diagrams 
also  are  to  be  interpreted  as  approximations  (in  this  and  the  following 
figures,  the  value  3i  which  is,  according  to  equation  (l),  physically 

more  important  and  which  characterizes  the  logarithmic  increment  of  the 
excited  perturbation  amplitude  has  been  plotted  instead  of  c±) .  Here, 
as  in  Schlichting's  report  (ref.  k)   mentioned  before,  only  the  derivative 


Pi  8* 


d(cc&*)  \  U 
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at  the  location  =  0,  that  is,  the  slope  of  the  -rr —  curve  at 

TT  a 

ua 

the  ccS*-axis  is  rigorously  correct.   H.  Schlichting  interpolated 

between  the  base  points  at  which  he  had  determined  the  slope  directly, 

although  with  a  much  higher  calculation  expenditure,  with  a  curve  of  the 

third  degree 

— —  =  aQ  +  a1(aS*)  +  a2(aS*)2  +  a^aS*)3 
ua 

where  the  four  constants  are  fixed  by  the  coordinates  of  the  base 
points  and  the  values  of  the  curve  slope  in  them.   Thereby, 

Pi  5* 
Schlichting  obtained  for  the  Blasius  profile  higher  values  — - —  than 

ua 

those  occurring  in  figure  25  whereas  the  values  of  the  slope  are  in 

3-5* 
agreement.   Actually,  the  curves  — - — (ct5*)  in  the  center  part  of  the 

Ua 

(aS*) -region  enclosed  in  each  case  will  run  somewhat  higher  or  lower 

than  indicated  by  figures  22  to  26.   At  any  rate,  however,  they  may  be 

interpreted  as  a  first  approximation  in  the  usual  sense. 

3iS*         * 
These  curves  — —  over  ct&   represent  sections  through  the 

Ua 

"excitation  mountain  range"  enclosed  by  the  neutral  curve  as  base  curve, 

along  the  lines  yc  =  Const,   or  cr  =  Const.   By  interpolation  one 

obtains  from  them  the  maps  of  excitation  represented  in  figures  27  to  32, 

3n-S* 
in  which  the  lines  of  equal  excitation  --= —  =  Const.,  can  be  interpreted 

Ua 
as  "contour  lines"  of  the  "excitation  mountain  range."  Instead  of  the 
lines  of  intersection  cr  =  Const.,  we  plotted  in  figures  27  to  32  the 
lines 

p  v   crao* 

—r-   =  —  =   Const. 

Ua2   UaRe* 

the  significance  of  which  will  be  discussed  later  in  section  VII.   Even 
at  first  consideration  of  these  excitation  maps,  a  fundamental  difference 
in  the  shape  of  the  excitation  mountain  range  is  conspicuous  according 
to  whether  the  velocity  profiles  of  the  laminar  friction  layer  lie  in 
the  region  of  decreasing  pressure  (3  >  0)  or  of  increasing  pressure 
(3  <  0). 

In  the  region  of  decreasing  pressure,  the  "excitation  mountain 
range"  has  the  form  of  a  mountain  with  pronounced  peak  which  is  steeply 
ascending  for  a  small  Reynolds  number  Re*,  slowly  flattens  after  a 
larger  Re*  and  shrinks  to  zero  width  and  height  for  Re* — ><».   The 
absolute  height  of  the  peak,  that  is,  the  maximum  excitation  increases 
with  decreasing  3- 
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In  the  region  of  increasing  pressure,  in  contrast,  the  excitation 
mountain  range  changes  behind  the  peak  with  growing  Re*   into  a 
"mountain  ridge"  of  constant  width  and  constant  contour  profile. 

The  properties  of  this  "ridge"  have  been  thoroughly  investigated 
by  W.  Tollmien  (ref.  5)-   Since  we  made  use  of  Tollmien's  theory,  we 
must  now  briefly  represent  its  results. 

Searching  for  a  general  instability  criterion,  W.  Tollmien  estab- 
lished that  the  frictionless  perturbation  equation  (5)  possesses  for  the 
laminar  velocity  profiles  in  the  region  of  pressure  increase  -  which 
have  a  point  of  inflection  in  contrast  to  those  of  the  region  of 
decreasing  pressure  -  for  Re* — ►  °°  aside  from  the  neutral  solution 
existing  for  all  profiles  with  the  parameters 


a  =  0 


<PnI 


U 


(93) 


in  addition,  the  neutral  solution  with  the  parameters 


a  =  etc 


cr  =  U£ 


fPnII  =  9s 


(9*0 


the  subscript  s  therein  denotes  the  point  of  inflection. 

W.  Tollmien  calculated  the  eigenvalue  as  and  the  eigenfunction 
qps  for  the  sinusoidal  velocity  distribution  (30)  which  we  also  used 
from  the  frictionless  perturbation  equation  (5) 


(aa). 


p  cot  p 


sin  p   2(1   -   £ 


(95) 


9S  = 


y 

sin  p  — 

^  a 


sin 


?*  = 


ps 


sin  p   aafl  -  — 

e   \    a 


sin 


ps 


(0  <  y  <  a) 


(y>a) 


(96) 


3^ 
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For  this  second  neutral  eigensolution,  the  phase  velocity  cr  therefore 
equals  the  velocity  of  the  basic  flow  at  the  point  of  inflection  Us. 

For  the  excitation  $±     and  the  circular  frequency  3r   in  the 
neighborhood  of  the  neutral  frequencies  a&*  =  0  and  aS*  =  (aS*)s, 
W.  Tollmien  has  derived  the  following  formulas 


where  the  subscript  w  signifies  that  the  values  have  to  be  taken  at 
the  wall  y  =  0,  and  that 


a  =  as:      — —  =  B 


[aS*)3    -    (as&*)2a5*|,    ^  Re 


Us      R*        BE  Us 

=  ul  a&    "  ~  uT7^" 


(a5*)3  -    (as5*)2ao*j 


(98) 


with 


where 


B  =  n 


lim 

€  — »0 


Uc 


Uc 


<Psdy 


ua& 


*2 


U 


0        (U   -  Us) 


E^^L 


9o2dy  + 


2  ^s 


Uc 


+e    (U   -  Us)2  "S 


s+e 


cp  2dy 


(99) 


(100) 


We  have  calculated  the  "ridge  contour  profiles"  of  the  excitation 
and  the  circular  frequency  according  to  the  formulas  (97)  and  (98)  for 
the  Hartree  velocity  distributions  3  =  -0.10,  -O.lU,  -0.l6,  -0.18, 
-O.I98;  the  values  for  the  derivatives  of  the  velocity  were  not  taken 
from  the  sinusoidal  approximation  distribution  but  from  the  exact 
velocity  distribution.   These  "ridge  contour  profiles"  of  the  excitation 
and  the  circular  frequency  have  been  plotted  in  figures  33  and  3k.      The 
curves  corresponding  to  the  formulas  just  mentioned  have  been  drawn  in 
solid  lines j  the  transitions  between  the  two  curve  arcs,  interpolated 
somewhat  arbitrarily,  have  been  drawn  in  dashed  lines.   In  the  range 
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of  very  large  Re"*,  the  lines  of  constant  excitation  and  constant  cir- 
cular frequency  were  transferred  from  figures  33  and  3^  into  figures  31 
and  32. 

Figure  35  represents  the  variation  of  the  maximum  excitation 

\  Ua  /max 

against  the  form  parameter  3>  separately  for  the  excitation  for  finite 
Re*   (peak  of  the  excitation  mountain  range)  and  for  Re*  >  °o  (ridge 
of  the  excitation  mountain  range).   For  3  <  -0.10  these  two  values 
seem  more  and  more  to  approach  one  another  so  that  with  decreasing  3 
the  peak  becomes  less  and  less  pronounced,  and  the  excitation  starting 
from  small  Re*  monotonically  increases  to  the  values  of  the  "ridge 
contour  profile." 

This  figure  and  the  preceding  ones  show  clearly  that  the  maximum 
excitation  and  accordingly  the  excitation  in  general  is  considerably 
larger  in  the  region  of  pressure  increase  than  in  that  of  pressure 
decrease  for  smaller  Re*  as  well. 

After  thus  having  estimated  the  magnitude  of  excitation  in  the 
entire  instability  range  of  all  velocity  distributions  occurring  in  the 
laminar  friction  layer,  we  shall  discuss  the  physical  conclusions 
resulting  from  our  calculations  for  the  position  of  the  transition  point. 

VII.  DISCUSSION  OF  THE  RESULTS 

Let  the  pressure  distribution  Ua(x)   against  the  arc  length  of 
the  cross  section  profile,  and  the  oscillation  of  circular  time  frequency 
3r  be  prescribed.   Then  this  oscillation  superposed  on  the  boundary- 
layer  flow  travels  downstream  on  a  curve 

3rv 

■^P  =  f(x) 

Ua2 

As  we  mentioned  at  the  beginning,  there  pertains  to  every  point  x  of 
the  profile  in  the  flow  a  fixed  value  of  the  Pohlhausen  parameter  X 
which  characterizes  by  way  of  approximation  the  velocity  distribution 
in  the  boundary  layer  at  this  point.  According  to  figure  3^  i-n  refer- 
ence 3  one  may  coordinate  to  this  parameter  X  the  Hartree  parameter  3 
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and  therewith  one  of  the  excitation  maps  calculated  in  the  present 
report  (figs.  27  to  32) .  In  order  to  obtain  the  excitation  of  the 
perturbation  3r  at  a  certain  point  x  of  the  profile,  one  has  there- 

3-&* 


fore  to  read  off  the  excitation 


Uc 


'number  of  the  "contour  line") 


on  the  corresponding  excitation  map  at  the  point  of  the  map  determined 
by  the  pair  of  values  of  the  Reynolds  number 


Re 


*   UaS* 


;(x) 


and  the  dimensionless  circular  frequency 


3rv 


=  f  x 


Us 


We  shall  start  the  discussion  of  the  results  of  our  above  excitation 

Uoot 

calculations  with  the  limiting  case  of  a  small  Reynolds  number  Re  =  . 

According  to  the  explanations  in  section  I  (compare  fig.  1)  the  stability 
limit  lies,  for  small  Reynolds  numbers  only,  at  the  point  where  the 
laminar  boundary  layer  separates.   If  the  perturbation  waves  are  so 
long  that  the  curve 


v 

U  2 


f(x) 


intersects  the  instability  region  for  the  separation  profile  (3  =  -0.198), 
it  is  very  violently  excited  when  entering  this  zone  and  leads  quickly 
to  transition  to  the  turbulent  flow  pattern.   The  transition  point  then 
practically  coincides  with  the  separation  point  which  we  had  denoted 
its  rearward  limit.   If,  on  the  other  hand,  the  perturbation  waves  are 
very  short  so  that  the  curve 


3rv 


f(x) 


does  not  intersect  the  instability  region,  the  laminar  layer  separates 
without  transition. 
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If  the  Reynolds  number  Re  =  — ^-  increases,  the  stability  limit 

shifts  forward  in  the  direction  toward  the  pressure  minimum  (fig.  1), 
the  perturbation  enters  a  region  of  instability  further  upstream  for  a 
larger  value  of  the  parameter  (3  or  the  Pohlhausen  parameter   A.,  and 
is  there  initially  excited  to  a  degree  which  decreases  the  more  the 
stability  limit  shifts  forward  but  which  increases  due  to  the  fact  that 
the  perturbation  downstream  (with  decreasing  p   or   A.)  reaches  insta- 
bility zones  with  rapidly  increasing  excitation.   The  transition  point 
shifts  frontward  corresponding  to  the  excitation  which  started  earlier 
and  is  still  strong. 

If  we  finally  increase  the  Reynolds  number  Re  =  -^—     so  that  the 

stability  limit  shifts  ahead  of  the  pressure  minimum,  the  excitation  in 
turn  starts  accordingly  sooner;  however  -  and  this  must  be  regarded  as 
the  most  important  result  of  our  calculations  for  the  time  being  -  in 
the  region  of  decreasing  pressure,  the  excitation  is  so  slight  that  it 
generally  attains  amounts  equalling  the  excitations  produced  in  the 
cases  treated  just  now  only  after  having  passed  the  pressure  minimum. 

In  this  manner,  one  may  easily  give  the  theoretical  explanation 
for  the  fact  proved  by  many  experiments,  that  the  transition  point  even 
in  case  of  very  high  Reynolds  numbers  rarely  ever  shifts  ahead  of  the 
pressure  minimum.   Only  in  cases  of  very  long  acceleration  sections 
and  perturbations  with  very  long  waves  where  the  perturbation  does  not 
too  soon  leave  the  instability  region  again,  the  limited  excitation  of 
the  region  of  decreasing  pressure  will  be  sufficient  to  induce  the 
transition  still  in  the  region  of  pressure  decrease.   However,  these 
cases  are  rare  in  technical  application. 

A  detailed  calculation  of  the  degree  of  excitation  which  causes 
the  transition  is  meaningful  only  in  connection  with  corresponding 
experiments.   Therefore,  it  will  be  postponed  until  these  experiments, 
now  in  the  preparatory  stage,  have  been  carried  out.   Probably  one  will 
have  to  regard  the  form  of  the  pressure  distribution  in  the  region  of 
pressure  increase  as  the  most  important  test  condition  since,  according 
to  our  theoretical  deliberations,  the  contributions  to  the  excitation 
of  perturbations  in  the  region  of  pressure  decrease  are  insignificant. 
The  aim  of  this  experimental  investigation  and  of  the  excitation  cal- 
culation to  be  performed  simultaneously  on  the  basis  of  the  present 
report  will  be  to  find  a  connection  between  the  pressure-distribution 
form  and  the  degree  of  excitation  attained  at  the  measured  transition 
point  so  that  it  will  be  possible  to  calculate,  inversely,  the  transition 
point  for  a  prescribed  pressure  distribution  from  this  relation. 
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VIII.  SUMMARY 


As  a  contribution  to  the  solution  of  the  important  problem  of  the 
calculation  of  the  transition  point  of  a  plane  laminar  flow,  we  had 
first  determined  (in  an  earlier  report,  according  to  Tollmien's  method 
of  small  oscillations  for  the  Hartree  velocity  distributions  appearing 
in  the  boundary  layer  in  case  of  decreasing  and  of  increasing  pressure) 
only  the  critical  Reynolds  number  beyond  which  the  perturbations  super- 
posed on  the  laminar  flow  are  excited.   In  connection  with  those  cal- 
culations now,  the  excitation  itself  in  the  entire  instability  range 
of  the  perturbations  was  calculated.   The  excitation  in  the  narrow 
instability  range  of  decreasing  pressure  turns  out  to  be  very  much 
smaller  than  the  excitation  in  the  more  extensive  instability  range  of 
increasing  pressure;  thus  the  known  fact  that  the  transition  point 
generally  does  not  shift  ahead  of  the  pressure  minimum  even  in  case  of 
high  Reynolds  numbers  may  be  explained  on  a  theoretical  basis,  as  shown 
in  tables  1,  2,  and  3- 

Systematic  experimental  measurements  of  the  transition  point, 
together  with  calculations  to  be  performed  on  the  basis  of  the  results 
given  here,  are  to  establish  the  connection  between  the  variation  of 
the  pressure  gradient  and  the  degree  of  excitation  which  produces  the 
transition  and  thereby  a  basis  for  determination  of  the  transition  point 
for  prescribed  pressure  variation  by  calculation. 


Translated  by  Mary  L.  Mahler 
National  Advisory  Committee 
for  Aeronautics 
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TABLE   1 
SERIES  COEFFICIENTS  OF  THE  SOLUTIONS     <Pj*,      tp  *     FOR  THE  BASIC  FLOW 
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TABLE   1 
SEBIES  COEFFICIENTS  OF  THE  SOLUTIONS     (p1#,      rp^'      FOB  THE   BASIC  FLOW 
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TABLE   1 
SERIES  COEFFICIENTS  OF  THE  SOLUTIONS     <Pj*,      <j>2*      FOR  THE  BASIC  FLOW 
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TABLE  2 
SERIES  COEFFICIENTS  OF  THE  SOLUTIONS     cp-,*,      <P2*     F0R  TEE  BASIC  FL0W 

U   =  0.    +   (1   -  U8)    sin  f^Ll-8.  2).      ABBREVIATION:      p  = 
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SERIES  COEFFICIENTS  OF  TEE  SOLUTIONS  tjlj* ,   <pg*  FOR  THE  BASIC  FLOW 
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TABLE  3 

VALUES  OF  THE  DIFFERENTIAL  QUOTIENTS  OF  THE  REAL  AND  OF  THE  IMAGINARY 

PART  OF  THE  FUNCTION  ^U    \     WITH  RESPECT  TO  t\q 


c3Fr 

*1 

5Fr 

uFi. 

^0 

St>0 

Slio 

^0 

^0 

Stio 

-2 

0.135 

-0.226 

-3-6 

0.396 

-0.066 

-2.1 

.121 

-.230 

-3-7 

•  399 

-.030 

-2.2 

.117 

-.235 

-3.8 

•  395 

-.010 

-2.3 

.118 

-.239 

-3-9 

•  379 

-.055 

-2.4 

.121 

-.242 

-4.0 

.360 

-.100 

-2.5 

.128 

-.245 

-4.1 

.3^0 

-.155 

-2.6 

.135 

-.246 

-4.2 

.297 

-.205 

-2.7 

.145 

-.246 

-4.3 

.250 

-.270 

-2.8 

•157 

-.243 

-4.4 

.180 

-.338 

-2.9 

.172 

-.237 

-4.5 

.080 

-.341 

-3.0 

.191 

-.226 

-4.6 

.020 

-.282 

-3.1 

.219 

-.210 

-4.7 

-.026 

-.226 

-3.2 

.256 

-.188 

-4.8 

-.045 

-.182 

-3.3 

.296 

-.162 

-4.9 

-.057 

-.150 

-3-4 

.340 

-.133 

-5-0 

-.062 

-.118 

-3-5 

.380 

-.102 
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Figures  2,  3,  and  4.-    Polar  diagrams  for  determination  of  the  curves  of  con- 
stant excitation. 
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Figures  5,6,  and  7.-    Polar  diagrams  for  determination  of  the  curves  of  con- 
stant excitation. 
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Figures  8,9,  and  10.-    Polar  diagrams  for  determination  of  the  curves  of  con- 
stant excitation. 
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Figures  11,  12,  and  13.-    Polar  diagrams  for  determination  of  the  curves  of 

constant  excitation. 
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Figures  14,  15,  and  16.-    Polar  diagrams  for  determination  of  the  curves  of 

constant  excitation. 
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Figures  17  and  18.- 


Polar  diagrams  for  determination  of  the  curves  of  con- 
stant excitation. 
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Figures  19,  20,  and  21.-    Polar  diagrams  for  determination  of  the  curves  of 

constant  excitation. 
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Figures  22,  23,  and  24.-    The  excitation    £i —    as  a  function  of  the  reciprocal 

Ua 

perturbation  wave  length    a  =  —    for  constant  critical  velocity    c  . 
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Figures  27  and  28.-    The  curves  of  constant  excitation  in  the  instability  regions 
of  a  few  laminar  boundary-layer  profiles. 
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Figure  32.-    The  curves  of  constant  excitation  in  the  instability  regions  of  a 
few  laminar  boundary -layer  profiles. 
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Figure  34.-    The  circular  time  frequency    — -  Re*    for  boundary -layer  pro- 
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